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A partial wave analysis of the centrally produced 7r°7r° channel has been performed in pp 
collisions using an incident beam momentum of 450 GeV/c. An unambiguous physical solution 



has been found. Evidence is found for the / (980), / (1300) and / (1500) in the the S-wave. 
and the / 2 (1270) is observed dominantly in the .D^-wave. In addition, there is evidence for a 
broad enhancement in the D-wave below 1 GeV. 
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The study of the 7r°7r° system produced in pp annihilations [[TJ and n induced interactions 0, 
|3| has revealed that the mv S-wave has a complicated structure which indicates the presence 
of several scalar states. Some of these states may have a significant gluonic component. In 
order to understand more fully the gluonic nature of these scalar mesons a systematic study is 
required using as many decay modes and production processes as possible. This paper presents 
a study of the 7r°7r° final state formed in central pp interactions which are predicted to be a 
source of gluonic final states via double Pomeron exchange |4| and gives new information which 
may help in understanding the scalar meson spectrum. 

The reaction 

pp -> p f (n°iT )p s (1) 

has been studied at 450 GeV/c. The subscripts / and s indicate the fastest and slowest particles 
in the laboratory respectively. The WA102 experiment has been performed using the CERN 
Omega Spectrometer, the layout of which is described in ref. ||. Reaction (1) has been isolated 
from the sample of events having two outgoing charged tracks and four 7s reconstructed in the 
GAMS-4000 calorimeter, by first imposing the following cuts on the components of the missing 
momentum: |missing P x \ < 10.0 GeV/c, |missing P y \ < 0.25 GeV/c and |missing P z \ < 0.25 
GeV/c, where the x axis is along the beam direction. A correlation between pulse-height and 
momentum obtained from a system of scintillation counters was used to ensure that the slow 
particle was a proton. 

Fig. 0a) shows the two photon mass spectrum for 47-events when the mass of the other 
7-pair lies within a band around the 7r° mass (85-185 MeV). A clear rc° signal is observed with 
a small background. Events belonging to reaction (|l|) have been selected using a kinematical 
fit (6C fit, four-momentum conservation being used and the masses of two 7r°s being fixed). By 
requiring that the x 2 value for the 7r°7r° hypothesis is x 2 < 12.6, a total of 208 452 events have 
been selected. 

The J0/7T effective mass spectrum (fig. |l]b)) shows a clear A + (1232) signal which has been 
removed by requiring M(pfn°) > 1.5 GeV. There is no evidence for A production in the p s 7T° 
effective mass spectrum (not shown). 

The resulting centrally produced 7r°7r° effective mass distribution is shown in fig. [TJc) and 
consists of 166 597 events. A peak corresponding to the ^(1270) and a sharp drop at 1 GeV 
can be observed. 

A Partial Wave Analysis (PWA) of the centrally produced 7r°7r° system has been performed 
assuming the 7r°7r° system is produced by the collision of two particles (referred to as exchanged 
particles) emitted by the scattered protons. The z axis is defined by the momentum vector 
of the exchanged particle with the greatest four-momentum transferred in the centre of 
mass. The y axis is defined by the cross product of the two exchanged particles in the pp centre 
of mass. The two variables needed to specify the decay process were taken as the polar and 
azimuthal angles (9, 0) of the ir° in the 7r°7r° centre of mass relative to the coordinate system 
described above. 

The acceptance corrected moments \f^nt LM , defined by 

J(ft) = $> y L °(fi) + 2 £ t LMJ Re{F L M (fi)} (2) 

L L,M>0 
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have been rescaled to the total number of observed events and are shown in fig. ^|. The moments 
with M > 2 (i.e. £43 and £44) and all the moments with L > 4 (not shown) are small and hence 
only partial waves with spin I = and 2 and absolute values of spin ^-projection m < 1 have 
been included in the PWA. 

An interesting feature of the moments is the presence of some structure in the L = 2 and 
L = 4 moments for 7r°7r° masses below 1 GeV which indicates the presence of D-waves. This 
type of structure has not been observed in rx induced reactions [§, ||. In order to see if this 
effect is due to acceptance problems or problems due to non-central events, we have reanalysed 
the data using a series of different cuts. Firstly, we required that M(pfTX°) > 2.0 GeV; however, 
it was found that after acceptance correction the moments were compatible with the set for 
M(pf7r°) > 1.5 GeV, showing that diffractive resonances in the range 1.5 < M(p f Tx°) < 2.0 GeV 
have a negligible effect on the moments. We have also required that the rapidity gap between 
any proton and tx° in the event is greater than 2 units. Again the resulting acceptance corrected 
moments did not change. In order to investigate any systematic effects we have also analysed 
the central tx + tx~ data and a similar structure is also found in the L = 2 and L = 4 moments |J. 

It is interesting to compare the above result with previously published data on other centrally 
produced txtx systems. In preliminary results from the E690 experiment at Fermilab, activity 
at a similar level is observed in the £40 moment below 1 GeV ]7). There is also evidence for 
this structure in the 7r°7r° data from the CERN NA12/2 experiment ||. The AFS experiment 
at the CERN ISR also observed moments that the moment deviated from zero in this mass 
region; however, in their analysis they claimed this deviation was due to problems of the Monte 
Carlo simulating low energy tracks ||. 

This structure does indeed seem to be a real effect which is present in centrally produced 
txtx systems. It has recently been suggested |H| |TTJ that central production may be due to the 
fusion of two vector particles and this may explain why higher angular momentum systems can 
be produced at lower masses. 



The amplitudes used for the PWA are defined in the reflectivity basis |12|] . In this basis the 



angular distribution is given by a sum of two non-interfering terms corresponding to negative 
and positive values of reflectivity. The waves used were of the form with J = S and D, 
m = 0, 1 and reflectivity e = ±1. The expressions relating the moments (£lm) and the waves 
( J^J are given in table [I]. Since the overall phase for each reflectivity is indeterminate, one wave 
in each reflectivity can be set to be real (Sq and Df for example) and hence two phases can 
be set to zero {4> s - and cf) D + have been chosen). This results in 6 parameters to be determined 
from the fit to the angular distributions. 

The PWA has been performed independently in 40 MeV intervals of the 7r°7r° mass spectrum. 
In each mass an event-by-event maximum likelihood method has been used. The function 

N 

F = -J2 ln{I(Q)} + £ t LM e LM (3) 

i=l L,M 

has been minimised, where N is the number of events in a given mass bin, e^M are the effi- 
ciency corrections calculated in the centre of the bin and t^u are the moments of the angular 
distribution. The moments calculated from the partial amplitudes are shown superimposed on 
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the experimental moments in fig |2|. As can be seen the results of the fit well reproduce the 
experimental moments. 

The system of equations which express the moments via the partial wave amplitudes is 
non-linear which leads to inherent ambiguities. For a system with S and D waves there are 
two solutions for each mass bin. In each mass bin one of these solutions is found from the fit 
to the experimental angular distributions, the other one can then be calculated by the method 
described in ref. JL2]. In the case under study the bootstrapping procedure is trivial because 
the Barrelet function has only two roots and their real and imaginary parts do not cross zero 
as functions of mass, as seen in fig ||]a) and b). In order to link the solutions in adjacent mass 
bins, the roots are sorted over real parts in each bin in such a way that the real part of the first 
root should be larger than the real part of the second root (real parts of the two roots have 
different signs). For the first solution the imaginary parts of both roots are taken positive, the 
second solution is obtained by complex conjugation of one of the roots. 

The two PWA solutions are shown in fig. [|. By definition both solutions give identical 
moments and identical values of the likelihood. The only way to differentiate between the 
solutions, if different, is to apply some external physical test, such as requiring that at threshold 
the S-wave is the dominant wave. For one solution the D-waves are the dominant contribution 
near threshold. This solution has been rejected as unphysical. The S-wave for the physical 
solution is characterised by a broad enhancement below 1 GeV and two shoulders, at 1 and 1.4 
GeV. Broad enhancements are also seen in the three D-waves at low mass, but the intensities 
of the D-waves near threshold are several times smaller than that of the S-wave. A peak 
corresponding to the /2(1270) is clearly seen in the Dq wave, such a peak is less prominent in 
the Di wave and is absent in the Df wave. 

In order to obtain a satisfactory fit to the Sq wave from threshold to 2 GeV it has been 
found to be necessary to use three interfering Breit-Wigners to describe the / (980), / (1300) 
and /o(1500) and a background of the form a(m — m t h) b exp(—cm — dm 2 ), where m is the 
7T 7r° mass, m t h is the 7r°7r° threshold mass and a, b, c, d are fit parameters. The Breit-Wigners 
have been convoluted with a Gaussian to account for the experimental mass resolution. The 
fit is shown in fig. |3]c) for the entire mass range and in fig. |3]d) for masses above 1 GeV. 
The parameters used to describe the /o(1300) and /o(1500) are those found from a fit to the 
7r + 7r~ mass spectrum of the same experiment ||, namely 

/ (1300) M = 1308 MeV, T = 222 MeV 
/ (1500) M = 1502 MeV, T = 131 MeV. 

The parameters found for the /o(980) are 

/o(980) M = 986 ± 10 MeV, T = 76 ±15 MeV, 

and as can be seen, the fit well describes the 7r°7r° S^-wave spectrum. 

In conclusion, a partial wave analysis of a high statistics sample of centrally produced 
7T 7r° events has been performed. An unambiguous physical solution has been found. The S- 
wave is found to dominate the mass spectrum and is composed of a broad enhancement at 
threshold, a sharp drop at 1 GeV due to the interference between the / (980) and the S-wave 
background, the / (1300) and the / (1500). The D-wave shows evidence for the / 2 (1270) and 
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a broad enhancement below 1 GeV. It is interesting to note that the /2(1270) is produced 
dominantly with m = 0. 
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Table 1: The moments of the angular distribution expressed in terms of the partial waves. 



V47rt o = 


l^o \ +\D \ 2 + \D x | 2 + \D+\ 2 


V&rho = 


^(2|D ~| 2 + |Dr| 2 + \D+\ 2 ) + 2|S- -||Do | cos(0 5 - - <j> D -) 


V&rt 21 = 


^PrlPol cos(0 D - - <f> D -) + v^i^lPrl cos(0 s<7 - <f, D -) 




2 m\ 2 -\Dt\>) 


VAttUo = 


f|D -| 2 -KIDrl 2 + l^ + l 2 ) 


V / 47Tt 4 l = 


^polPrloos^--^-) 


V / 47T^42 = 


^m\ 2 -\Dt\ 2 ) 
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Figures 

Figure 1: a) Effective 77 mass for 47-events when the mass of the other 7-pair lies within the 
7T° mass band, b) the M(p/7r°) and c) the centrally produced 7r°7r° effective mass spectrum. 



Figure 2: The \f^nt LM moments from the data. Superimposed as a solid histogram are the 
resulting moments calculated from the PWA of the 7r°7r° final state. 



Figure 3: a) The real and b) Imaginary parts of the roots (see text) as a function of mass 
obtained from the PWA. c) and d) The 7r°7r° Sq wave with fit described in the text. 



Figure 4: The physical (solid circles) and unphysical (open circles) solutions from the PWA 
of the 7r°7r° final state. 
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Figure 1 
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Figure 3 
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